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Abstract Ttie increasing level of experimental control over 
atomic and optical systems gained in the past years have paved 
the way for the exploration of new physical regimes in quantum 
optics and atomic physics, characterised by the appearance of 
quantum many-body phenomena, originally encountered only 
in condensed-matter physics, and the possibility of experimen- 
tally accessing them in a more controlled manner In this review 
article we survey recent theoretical studies concerning the use 
of cavity quantum electrodynamics to create quantum many- 
body systems. Based on recent experimental progress in the 
fabrication of arrays of interacting micro-cavities and on their 
coupling to atomic-like structures in several different physical 
architectures, we review proposals on the realisation of paradig- 
matic many-body models in such systems, such as the Bose- 
Hubbard and the anisotropic Heisenberg models. Such arrays 
of coupled cavities offer interesting properties as simulators of 
quantum many-body physics, including the full addressability 
of individual sites and the accessibility of inhomogeneous mod- 
els. 




Crystal light: Photons (yellow) that ordinarily move freely 
through an array of cavities (top) might be frozen in place by 
their mutual repulsion (bottom). Such repulsion can be induced 
by laser pumping of suitable atoms (red) placed in the cavities. 



> 

(N 
00 

o 

oo 

o 



Copyright line will be provided by the publisher 



Quantum Many-Body Phenomena in Coupled Cavity 
Arrays 

Michael J. Hartmann ^-^^s,*^ Fernando G.S.L. Brandao and Martin B. Plenio ^'^ 

Technische Universitat MUnchen, Physik Department I,James Franck Str., 85748 Garching, Germany 
J Institute for Mathematical Sciences, Imperial College London, SW7 2PG, United Kingdom 
QOLS, The Blackett Laboratory, Imperial College London, Prince Consort Road, SW7 2BW, United Kingdom 

Received: . . . 
Published online: . . . 



X 



'Key words: quantum many-body models, polaritons, cavity QED, photon blockade 
PACS: 42.50.Dv,42.50.Pq,03.67.Bg,75.10.Jm 



Contents 



10. Summary |24 



1 . Introduction 

2. Cavity QED 

3. Coupled Cavity Arrays 

4. The Bose-Hubbard Model in Coupled Cavity Arrays 

5. Models Based on the Jaynes-Cummings Hamiltonian 

6. Phase Diagrams 

7. Effective Spin Models in Coupled Cavity Arrays . . 

8. Promising Platforms for Realisations 

9. Applications 



[3 
3 

i 

13 
13 



1. Introduction 

Quantum optics and atomic physics are concerned with the 
interaction of light and matter in the quantum regime, typ- 
ically on the scale of a few atoms and light quanta. As 
atoms do not naturally interact strongly with each other 
and with light, the physics of such systems can usually 
be very well understood by ignoring collective effects and 
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treating the interactions perturbatively. This is a much sim- 
pler and cleaner situation than the one encountered in the 
condensed matter context. There, strong interactions among 
the basic constituents, such as nuclei and electrons, lead to 
the emergence of completely new physics when one con- 
siders a mesoscopic or macroscopic number of interacting 
particles. Thus, even though the fundamental interactions 
between the constituent particles are usually known, it is 
challenging to fully describe the properties of such sys- 
tems. 

The higher level of isolation of quantum optical sys- 
tems have permitted, over the past decades, the develop- 
ment of experimental techniques for their manipulation 
with an unprecedented level of control. Highlights include 
the invention of trapping and cooling techniques [T], which 
have found applications e.g. in spectroscopy and precision 
measurements and which culminated in the realisation of a 
Bose-Einstein condensate (BEC) of neutral atoms |2l; the 
manipulation of the state of individual atomic and photonic 
systems at the quantum level IMM and its use in the study 
of the foundations of quantum theory |4l!5l; and more re- 
cently the realisation of primitives for the processing of 
quantum information and for the implementation of quan- 
tum computation in several different quantum atomic and 
optical set-ups IS). 

This progress raised the possibility of engineering strong 
interactions in these systems and of realising such strongly 
correlated models in large scales. This is of interest from a 
fundamental point of view, as in this way we can form new 
phases of matter for atoms, photons, or even molecules 
that do not exist outside the laboratory. Moreover, due to 
the high level of experimental control over them, it is pos- 
sible to study many-body physics in a much cleaner man- 
ner than it is possible employing condensed-matter sys- 
tems, in which the strong interaction with the environment 
and the very short space- and time-scales involved inhibit 
a precise control over their static and dynamic properties. 
The use of one quantum mechanical system to simulate 
the physics of a different one is an idea that dates back 
to Feynman Q and has attracted considerable interest re- 
cently in the field of quantum information science 161 . 

A very successful example in this direction is the study 
of cold trapped atoms in optical lattices |8]. By loading a 
BEC of neutral atoms in a periodic optical lattice, formed 
by the off -resonant dipole interaction of atoms with over- 
lapping standing waves of counter propagating lasers, the 
otherwise weak interactions of the atoms can be increased 
to such an extend that a strongly correlated many-body 
model is formed This idea led to several important 
experiments, such as the realisation of the Bose-Hubbard 
model 111 01 and of a Tonks-Girardeau gas ifTTI and is cur- 
rently an active area of theoretical and experimental re- 
search una. 

Also of note is the early study of strongly correlated 
many -body models in Josephson junctions arrays 1 13|. Al- 
though a Josephson junction is a mesoscopic solid state 
device, it behaves in many respects just like the type of 
systems studied in quantum optics, as the relevant physics 



can be understood by looking at a few quantum levels. 
In this context, arrays of interacting Josephson junctions 
have been used to reproduce properties of bosonic parti- 
cles. For example, the Mott insulator phase lfT4l and the 
superconductor-to-Mott insulator transition fl5\ have been 
observed in Josephson junctions arrays. 

Finally, theoretical proposals for using trapped ions for 
realising many-body models have been put forwarded 1161 
[TTlfTslfTQl and a first benchmark experiment in this direc- 
tion has been reported recently 1201 . 

This article is devoted to review theoretical proposals 
for the use of the atom-light interaction in coupled micro- 
cavity arrays to create strongly correlated many -body mod- 
els. This research was initiated in |21] and shortly after in 
122112311 (see also 1241 ). Since then a substantial number of 
publications have extended these ideas in a variety of di- 
rections. These include photonic Mott insulators II25II26I . 
spin models l27ll28l l29ll^0ll3T]l32l. dynamical effects l33l 
1 34, 35, 36, '37), multi component models ll38l . phase dia- 
grams 139, 40l !4ni42ll43l . alternative physical realisations 
l44B5ii46ii47il48l . experimental signatures F49IJ50J and quan- 
tum information applications |5 1 , 52 , 53 , 54 , 55 56j . In fact, 
historically, some quantum information proposals ISTl were 
presented before the theory for effective many-body mod- 
els EDEaiia. 

Recent experimental progress in the fabrication of micro- 
cavity arrays and the realisation of the quantum regime 
in the interaction of atomic-like structures and quantised 
electromagnetic modes inside those II57II5 8ll59ir60ll6 1 Il62l 
I63ll64l opened up the possibility of using them as quan- 
tum simulators of many -body physics. The first motivation 
for this study is the desire to make photons or combined 
photonic-atomic excitations enter new quantum states, which 
do not naturally appear in nature. The second is that this 
set-up offer advantages over other proposals for the reali- 
sation of strongly interacting many-body models in quan- 
tum optical systems. Due to the small separations between 
neighbouring sites, it is experimentally very challenging to 
access individual sites in Josephson junctions arrays and 
optical lattices. Neighbouring sites in coupled cavity ar- 
rays, on the other hand, are usually separated by dozens 
of micrometres and can therefore be accurately accessed 
by optical frequencies. This allows e.g. the measurement 
of local properties as well as the exploration of inhomoge- 
neous systems. 

As the field of coupled cavity arrays is still quite young 
and many preprints are currently being reviewed, we have 
restricted ourselves here to only present details of work 
that has been published in peer reviewed journals before 
completion of the present text. This has been done in order 
to avoid any interferences with the policies of journals that 
consider the work in question for publication. 

The organisation of the paper is as follows: In section 
|2]we discuss the use of cavities to enhance the interac- 
tion of light and matter and in section [3] we develop the 
basic description for an array of interacting cavities. The 
following section previews the theory for engineering ef- 
fective Bose-Hubbard Hamiltonians with polaritons and 
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Figure 1 Cavity QED: A two level atom interacts via a dipole 
coupling g with the photons in the cavity. Excitations a lost via 
spontaneous emission at a rate 7 and cavity decay at a rate k. 



photons. In section [5] we discuss a second class of effec- 
tive Hamiltonians that exhibits a phase diagram similar to 
the Bose-Hubbard model: Jaynes-Cummings based many 
body models. The phase diagrams of both classes are then 
reviewed in section |6] Effective spin lattice Hamiltonians 
can also be generated in coupled cavity arrays and we re- 
view these approaches in section [T] In section [8] we then 
discuss several experimental platforms which appear to be 
suitable for a future implementation of the theory approach 
presented in this review. We finish with a summary in sec- 
tion[T0l 



2. Cavity QED 

Usually photons interact only over a very short time, and 
hence very weakly, with atoms. This interaction can be 
substantially increased if the atoms and photons are trapped 
inside a resonator. In the simplest situation, we can con- 
sider the coupling of a single light mode with a single two 
level system inside a cavity. Using the dipole and the ro- 
tating wave approximations their dynamics is governed by 
the Jaynes-Cummings Hamiltonian, 

H = huJca^a + fi^o\e){e\ + hg{a^\g){e\+a\e){g\), (1) 

where ujc , ^0 the frequencies of the resonant mode of 
the cavity and of the atomic transition, respectively, g is 
Jaynes-Cummings coupling between the cavity mode and 
the two level system, is the creation operator of a photon 
in the resonant cavity mode, and | g) , |e) are the ground and 
excited states of the two level system. The Rabi frequency 
g is given by 

5 = d/(ro)./^I^, (2) 

where VInode is the mode volume of the cavity, /(rp) is the 
mode-function at the position tq of the two level system, 
d is the electric dipole transition moment, and eo is the 
free-space permittivity. 

In addition to the coherent interaction, there are two 
main loss processes that affect the dynamics of the sys- 
tem: spontaneous emission from level |e) to level | g) at a 



rate 7 and leaking out of photons of the cavity mode at rate 
K, see figure [U The cavity decay rate k is connected to the 
quality factor of the cavity Qhy Q = uJc / (2k). Although 
other processes, such as thermal motion of the atom or 
dephasing due to background electromagnetic fields, also 
contribute to the losses of the system, their effect are usu- 
ally much smaller and can be disregarded in many appli- 
cations. 

The strong coupling regime of cavity QED is reached 
when the cooperativity factor, given by the the vacuum 
Rabi frequency to the square over the product of the spon- 
taneous emission and cavity decay rates, 

^■■^1^, (3) 

is much larger than one: ^ ^ 1. In this regime the coher- 
ent part of the evolution dominates over the decoherence 
processes and quantum dynamics of the joint atom cav- 
ity mode system can be observed. Such a regime has first 
been achieved in seminal experiments for microwave ID 
and later optical frequencies 1651 . both using a single atom 
inside a Fabry-Perot cavity formed between two minia- 
ture spherical mirrors. In these works, the strong coupling 
regime has been employed e.g. to create conditional quan- 
tum dynamics and entanglement between the photonic and 
atomic degrees of freedom B1 l66l and to perform quantum 
non-demolition measurements (4\ . Recently a strong cou- 
pling regime with very high cooperativity factors has also 
been achieved with artificial atoms formed by supercon- 
ducting structures that interact with the electromagnetic 
modes of stripline resonator |62|. 

Most importantly for our purposes the strong coupling 
regime led to a realisation of the photon blockade effect 
[67 1, where, due to the strong interaction of the cavity 
mode with atoms, a single photon can modify the reso- 
nance frequency of the cavity mode in such a way that a 
second photon can not enter the cavity before the first leaks 
out. 

In this review we are interested in the situation where 
several cavities operating in the strong coupling regime 
are coupled to each other It turns out that the Fabry-Perot 
architecture is not very suitable for the coupling of sepa- 
rate cavities. Nonetheless, several new cavity QED set-ups 
have recently emerged in which (i) the strong coupling 
regime has been achieved and (ii) the construction of ar- 
rays of coupled cavities has been realised, or at least seem 
reasonable to be so. These new systems, which we briefly 
review in section[8] are the strongest candidates for realis- 
ing the proposals that we review in this paper 



3. Coupled Cavity Arrays 

In this section we present the theoretical description for 
the coupling of cavities in an array. Photon hopping can 
occur between neighbouring cavities due to the overlap of 
the spacial profile of the cavity modes, see figure |2] . In 
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Figure 2 An array of coupled cavities. Photon hopping occurs 
due to the overlap (shaded green) of the light modes (green lines) 
of adjacent cavities. 



order to model such a process, we follow references 11211 
l68ii69J and describe the array of cavities by a periodic di- 
electric constant, e(r) = e(r + Rn), where r is a given 
three dimensional vector, i? is a constant and n labels tu- 
pels of integers. In the Coulomb gauge the electromag- 
netic field is represented by a vector potential A satisfying 
V • (e(r) A) = 0. We can expand A in Wannier functions, 
tjR, each localised in a single cavity at location R = Rn. 
We describe this single cavity by the dielectric function 
eR(r). Then, from Maxwell's equations one can derive the 
eigenvalue equation 

^SM^WR - Vx (Vxwr)=0, (4) 

where the eigenvalue ujq is the square the resonance fre- 
quency of the cavity, which is independent of R due to the 
periodicity. Assuming that the Wannier functions decay 
sufficiently fast outside the cavity, only Wannier modes of 
nearest neighbour cavities have a non-vanishing overlap. 

Introducing the creation and annihilation operators of 
the Wannier modes, and aR, the Hamiltonian of the 
field can be written as 

H = u;cY^ (a^aR + 

R, 

+ 2ujca (ak«R.' + h-c) . (5) 

<R,R'> 

Here^ 

<R R'> '^he sum of all pairs of cavities which are 
nearest neighbours. Since a <C 1, we neglected counter 
rotating terms which contain products of two creation or 
two annihilation operators of Wannier modes in deriving 
(|5]l. The coupling parameter a is given by ||69l , 

a = j d\ (eR(r) - e(r)) w^wr, ; |R - R'| = i? , 

(6) 

and has been obtained numerically for specific models ||70|. 

In deriving Hamiltonian Q we have assumed that all 
the cavities have the same resonant frequency and that the 
overlap of the Wannier modes a is the same for all cavity- 
cavity interactions. In practise of course there will always 
be some disorder in the array and the resonance frequen- 
cies, wc'(R), and tunnelling rates, 2tJc(R)a((R) will dif- 
fer from cavity to cavity. The disorder in the array can even 
give rise to interesting effects such as the emergence of 
glassy phases, see|6] 



4. The Bose-Hubbard Model in Coupled 
Cavity Arrays 

Although realising Hamiltonian Q might be interesting 
for quantum information propagation |71|721, it is not the 
type of model that one would like to quantum simulate. In- 
deed, the model is harmonic and can be very easily solved 
in terms of Bloch waves, which allows for a simple under- 
standing of all its basic properties. The situation changes 
dramatically if we add an on-site interaction term. The in- 
terplay of tunnelling and interaction leads to interesting 
many-body physics. A paradigmatic model for such a sit- 
uation is the Bose-Hubbard Hamiltonian, which reads 

k <k,k'> 

+ Uj2bih{blbk-l), (7) 

k 

where &^ creates a boson at site fc, J is the hopping rate, U 
the on-site interaction strength, and /i the chemical poten- 
tial. This model, which was first discussed by Fisher et al 
in Ref. 17311741 as the bosonic counterpart of the Hubbard 
model 1751 , contains two different phases at zero tempera- 
ture that are separated by a quantum phase transition where 
the ration U /J crosses a critical value. Let us briefly de- 
scribe the two phases qualitatively. 

When J ^ U, the tunnelling term dominates and the 
lowest energy state of the system is a condensate of de- 
localised bosons: the system is in a superfluid phase. If 
we start to increase a repulsive on-site interaction, adding 
more than one particle to one site requires energy that grows 
quadratically with the particle number. It will thus be en- 
ergetically favourable for the bosons to evenly distribute 
across all lattice sites and avoid moving around in order 
to minimise the repulsion. Above a specific value of U/J, 
the system ceases to be a superfluid and becomes a Mott 
insulator, with a well defined number of localised parti- 
cles in each siteQ. The superfluid-to-Mott insulator phase 
transition is of a quantum character and is driven by quan- 
tum fluctuations, as it occurs even at zero temperature. An 
order parameter for such a transition is the variance of the 
number of bosons in a given site, {{blbk — {bl.bk))'^)- In the 
superfluid phase, the number of bosons in each site fluctu- 
ates and thus the variance has a non-zero value. In the Mott 
insulator phase, in turn, the variance is close to zero, as the 
bosons tend to be localised in each site. 

As we saw in section |3] the coupling between neigh- 
bouring cavities naturally leads to the tunnelling term for 
photons in arrays of coupled cavities. Hence, for a reali- 
sation of a Bose-Hubbard model a way to effectively gen- 
erate an interaction term is needed. Indeed, up to harmonic 
terms which can easily be compensated for, the Bose-Hubbard 
interaction term has the form of a Kerr self non-linearity 

' In the case of incommensurate filling factors, i.e. if the num- 
ber of particles is not a multiple of the number of lattice sites, the 
system will always be in a superfluid phase. 
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Figure 3 An array of cavities as described by our model. Photon 
hopping occurs due to the overlap (shaded green) of the light 
modes (green lines) of adjacent cavities. Atoms in each cavity 
(brown), which are driven by external lasers (blue) give rise to an 
on site potential. Reprinted with permission from 1211 . 
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Figure 4 The level structure and the possible transitions of one 
atom, ujc is the frequency of the cavity mode, O is the Rabi fre- 
quency of the driving by the laser, gn and g24 are the parameters 
of the respective dipole couplings and 5, A and e are detunings. 
Reprinted with permission from 1211 . 



(^fc)^^fc- Although such an interaction naturally appears in 
some media as a result of a non-zero third order electric 
susceptibility, its effect is negligible on the level of indi- 
vidual quanta, which is a reason for the difficulty of real- 
ising nonUnear optics for individual quantum systems. Us- 
ing the enhanced light-matter interaction in cavity QED, it 
is possible to engineer much stronger nonlinearities. Intu- 
itively the strong interactions of the light mode with atoms 
inside the cavity, under particular circumstances, mediates 
strong nonlinear interaction among the photons of the cav- 
ity mode. The strength of the nonlinearity can be increased 
even further if instead of considering photons as the bosonic 
particles of the model, one considers polaritons, joint photonic 
atomic excitations. In the next section we review the pro- 
posal for realising the Bose-Hubbard Hamiltonian using 
such polaritons. 



4.1. Polaritonic Bose-Hubbard model 



The polaritonic Bose-Hubbard model is hosted in an array 
of cavities that are filled with atoms of a particular 4 level 
structure, which are driven with an external laser, see fig- 
ures|3]and|4] Thereby the laser drives the atoms in the same 
manner as in Electromagnetically Induced Transparency 
(EIT) ll76l : The transitions between levels 2 and 3 are cou- 
pled to the laser field and the transitions between levels 
2-4 and 1-3 couple via dipole moments to the cavity reso- 
nance mode. Levels 1 and 2 are assumed to be metastable 
and their spontaneous emission rates are thus negligible. 

It has been shown by Imamoglu and co-workers, that 
this atom cavity system can exhibit a very large nonlin- 
earity Il77ll78ll79l[80l . and a similar nonlinearity has been 
observed experimentally ||67l . 

Considering the level structure of figurelH in a rotating 
frame with respect to iJo = '''^ ^"'^^ 



the atoms in the cavity reads. 



UJc<J33 



ujc [a' a + i) + 
the Hamiltonian of 



N 



= J2 {'^'22 + 6<yis + + sHi) (8) 

N 



where cr^; = \kj){lj \ projects level I of atom j to level k 
of the same atom, ujc is the frequency of the cavity mode, 
i7 is the Rabi frequency of the driving by the laser and 1713 
and g24 are the parameters of the dipole coupling of the 
cavity mode to the respective atomic transitions. 

In a cavity array, N atoms in each cavity couple to the 
cavity mode via the interaction and photons tunnel be- 
tween neighbouring cavities as described in equation (|5]). 
Hence the full Hamiltonian describing this system reads 



1 

(9) 



R R ^ 

+ 2u!ca (^a^aR' + h.c 

<R,R'> 

and can emulate a Bose-Hubbard model for polaritons as 
we will see in the following. 

Assuming that all atoms interact in the same way with 
the cavity mode, the description can be restricted to Dicke 
type dressed states, in which the atomic excitations are de- 
localised among all the atoms. In the case where 1724 = 
and £ = , level 4 of the atoms decouples from the 
dressed-state excitation manifolds ||79l . The Hamiltonian 
([8j can then be expressed in terms of the following cre- 
ation (and annihilation) operators: 



pi 



= 4 (9SI2 - ^a.^ 



(10) 



A{A±S) 



QS- 



12 



ga 
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where g = VNgis, B = y/g^T^, A = V4B2 
'-'12 



a 21 and S- 



2_ v^JV 



The operators p\, and describe polaritons, quasi 
particles formed by combinations of atom and photon ex- 
citations. In the subspace spanned by symmetric Dicke 
states and in the limit of large atom numbers, ^ 1, 
they satisfy bosonic commutation relations, 



[P3.P1] = and 



Pj.Pi 



5ji for j, ; = 0, - 



(11) 



p\f,p\ andp!_ thus describe independent bosonic particles. 
In terms of these polaritons, the Hamiltonian dH) for 524 = 
and e = reads. 



324=0, e=0 



MopSpo + Ai+P+P++A*-P-P- , (12) 



where the frequencies are given by /ip = 0, /i+ = ((5 — 
A)/2and ^_ = {5 + A)/2. 

The polaritons only contain atomic contributions in 
the two metastable states 1 and 2 but not in level 3, which 
shows spontaneous emission. Hence, these polaritons do 
not populate radiating atomic levels and are thus called 
dark state polaritons ||76| . 

As we will discuss now, the dynamics of the dark state 
polaritons pq decouples from the dynamics of the remain- 
ing species p+ and p_ for a suitable parameter regime. 
In this regime the frequency ^0 is sufficiently separated 
from and such that neither the interaction with the 
atomic level 4 nor the photon tunnelling between cavities 
can induce any mixing between the polariton species po 
and p±. 



4.1.1. Polariton-polariton interactions and polariton 
tunnelling 

The polaritons p\ interact with each other via a self Kerr 
nonlinearity. To derive this interaction, one first needs to 
write the full Hamiltonian , dHJ, in the polariton basis, 
expressing the operators '^22 ^^'^ '^24 

terms of pj, and pL. The coupling of the polaritons to 
the level 4 of the atoms via the dipole moment 324 reads. 



524 



N 



h.c. 



gfi 
52 



-.g24^ [s\iPl 



h.c 



(13) 

where S\-^ ~ ^j=i '^ii- deriving ( fTSl ). a rotating 
wave approximation was used: In a frame rotating with re- 
spect to (fT2] |. the polaritonic creation operators rotate with 
the frequencies no, fJ.+ and Furthermore, the operator 
rotates at the frequency 2/io, (c.f. dD). Hence, pro- 
vided that 



|g24|, |e|, 1^1 < Ia^+I, Ia^-I 



(14) 



all terms that rotate at frequencies 2/io ^ (/^+ i M-) or 
fiQ ± fi+ or fiQ ± /i„ can be neglected, which eliminates all 
interactions that would couple pi and 5*14 to the remaining 
polariton species. 

For \g2A9i^ / B'^l <C the coupling to level 4 can be 
treated perturbatively. This results in an energy shift of 2U 
with 

524 

A 



(15) 



and in an occupation probability of the state of one 5I4 
excitation of — 2C/ /Z\, which determines an effective decay 
rate for the polariton p\ via spontaneous emission from 
level 4. Note that [/ > for Z\ < and vice versa. In a 
similar way, the two photon detuning e leads to an energy 
shift of e g^ for the polariton pj, which plays the role 
of a chemical potential in the effective Hamiltonian. 

Hence, provided ( fT4] i holds, the Hamiltonian for the 
dark state polariton p\ can be written as 



H,ff = U 



2 9 i 

(Po) +e^plpo, (16) 



in the rotating frame. 

Let us now look at an array of interacting cavities, each 
under the conditions discussed above. The first term of the 
Hamiltonian ^ has already been incorporated in the po- 
lariton analysis for one individual cavity. The second term, 
which describes photon tunnelling, transforms into the po- 
lariton picture via (fTOl i. To distinguish between the dark 
state polaritons in different cavities, we introduce the no- 
tation to label the polariton pj in the cavity at position 
R. The photon hopping translates to polariton hoppings 
according to. 



[2^ 



(17) 



+ "terms for other polariton species" . 

Contributions of different polaritons decouple due to the 
separation of their frequencies ijq, fi+ and As a con- 
sequence the Hamiltonian for the polaritons p^ takes on 
the form (jTji, with 



U as given by equation (fTSl l and /i = eg'^/B'^. 



(18) 



4.1.2. Polariton lifetime: spontaneous emission and cavity 
decay 

We now discuss the loss sources in the system. Levels 1 
and 2 are metastable and therefore have negligible decay 
rates on the relevant time scales. The decay mechanism for 
the dark state polariton then originates mainly from pho- 
tons that leak out of the cavity and from the small, but non- 
zero, population of level 4, which leads to spontaneous 



Copyright line will be provided by the publisher 











Laser & Photon. Rev. 1, No. 1 (2008) 


r 


||||[^|||||||^|i^evi«w3 1 


7 



emission. The resulting effective decay rate for the dark 
state polariton reads 



2)^2^74, 



(19) 



where k is the cavity decay rate, 74 the spontaneous emis- 
sion rate from level 4, n is the average number of pho- 
tons, and O the Heaviside step function. Assuming that 
5i3 = .924, the maximal achievable rate U /Fq can be read- 
ily seen to be 1713 / ^jAKO{'n — 2)74. Hence, for large co- 
operativity factors, ^ 3> 1, ratios of [//Iq ^ 1 can be 
achieved, which ensures that the dynamics of the Bose- 
Hubbard model can be observed under realistic conditions. 



4.1.3. The phase transition 

Of great interest is of course whether the quantum phase 
transition from a Mott insulator to a superfluid state could 
be observed with the present polariton approach. In a sys- 
tem with on average one polariton per cavity, the Mott in- 
sulator state is characterised by the fact that the local state 
in each cavity is a single polariton Fock state with van- 
ishing fluctuations of the polariton number The superfluid 
phase in contrast shows polariton number fluctuations. 

Figure |5] shows numerical simulations of the full dy- 
namics of an array of three cavities as described by equa- 
tion Q, including spontaneous emission and cavity decay, 
and studies the number of polaritons in one cavity/site, 
ni = {p\pi) for site I, and the number fluctuations, Fi = 

{{pjpi)'^) — (pjpi)'^. A comparison between the full model 
(|9|l and the effective model (|7]i is done by considering the 
differences in the occupation numbers, Sni = [«;](.avities ~ 
['^'Ibh' "iiid number fluctuations, SFi [^dcavities 
at each time step 1211 . Initially there is exactly one po- 
lariton in each cavity and parameters for toroidal micro- 
cavities from ISTl . (724 = 513 = 2.5 x lO^s— 1, 74 = 73 = 
1.6 X lO^s-i and k = 0.4 x lO^s"!, are assumed. As 
the system is driven from a Mott insulator state to a super- 
fluid state by ramping up the driving laser f2, the particle 
number fluctuations increase significantly. The numerics 
in figure |5] show good agreement of the full dynamics as 
described by equation ^ and the dynamics of the corre- 
sponding Bose-Hubbard model (|7]i and thus confirm the 
possibility of observing the Mott insulator-to-superfluid 
transition in such a system. The oscillatory behaviour is 
related to the fact that the initial state of the system is, due 
to the nonzero tunnelling rate J, not its ground state. This 
initial state has been chosen since it's preparation in an ex- 
periment is expected to be easier compared to other states 
and figure |5] thus shows the dynamics that is expected to 
be observed in an experiment. 

The phase diagram for this model and its dependence 
on the atom number N have been analysed in lf39ll . where 
also a glassy phase has been predicted for a system with 
disorder, see section |6] for details. 



A polariton model which generates an effective Lieb- 
Liniger Hamiltonian and is capable of reproducing the Tonks 
Girardeau regime was proposed in ||451 . The scheme em- 
ploys atoms of the same level structure as in figure |4] that 
couple to light modes of a tapered optical fibre. Further- 
more quantum phase transitions of Dicke-model in cavity 
QED have been analysed 1821 . 



4. 1 .4. Two component Bose-Hubbard model 

Before we conclude this section, let us briefly review the 
proposal of Ref. l38l for creating a two component polari- 
tonic Bose-Hubbard model. This model, which describes 
the dynamics of two independent bosonic particles in each 
site, reads 



H^ff - E '•'j 

k,j—b,c 

+ E 

kj—b,c 



<k.k' ^-.jjl—bjC 



1 +Et^^^«^(20) 



where bl,{cl) create polaritons of the type b{c) in the cav- 
ity at site k, n[,^'' = b^bk and n^'^'' = c\ck- ^^b and /ic are 
the polariton energies, Uh, Uc and Ub.c their on-site inter- 
actions and Jfc.b, Jc^c and J^.c their tunnelling rates. This 
model exhibits interesting quantum many-body phenom- 
ena such as spin density separation f83l, spin order in the 
Mott regime l84ll and phase separation |85|. 

The Hamiltonian ( l20l i can be realised in the same set- 
up as above, by considering the atom-photon interactions 
in a dispersive regime with 5 ^ ^, ff- Then, the polariton 
operators (doll read. 



P 



ga 



to leading order in S~^. Note that now, not only the po- 
lariton pq, but also pl_ do not experience loss from sponta- 
neous emission to leading order inS^^. One can therefore 
define two dark state polaritons species 



9^12 ^ 



ga' 



(21) 

From an analysis similar to the one outlined above for the 
one component case, it can be shown that the dynamics of 
these two polaritons species is given by the two component 
Bose-Hubbard model ( l20l i. where the parameters are given 



J - aK 

cc 



•he ~ OL^, Ub = 

2/2 r-,2N2 



2 2 
924g_i£_ 



by, Jbb 

2 2 r,2 2/2 r-,2N2 

TT — g24g ^{ a„H TL — 924 (ff ) p r- 

ther details of the derivation the reader is referred to Ref. 
f38l. 
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Figure 5 The Mott insulator to superfluid transition for 3 polaritons in 3 cavities compared to 3 particles in a 3 site Bose-Hubbard 
model, a: Log-plot of U and J and linear plot of the time dependent O (inset), b: Number fluctuations for polaritons in cavity I, Fi, for 
a single quantum jump trajectory, c: Difference between number fluctuations for polaritons in a cavity and number fluctuations in a pure 
Bose-Hubbard model, SFi, for a single quantum jump trajectory, d: Expectation value and fluctuations for the number of polaritons 
in one cavity/site according to the effective model Q with damping l |19l l. g24 = gis = 2.5 x lO^s— 1, 74 = 73 = 1.6 x lO^s"^, 
K = 0.4 X 10'^s"\ N = 1000, A = -2.0 x 10"'s"^ and 2uJca = 1.1 x 10'^s~\ The Rabi frequency of the driving laser is 
increased from initially J7 = 7.9 x 10^"s~^ to finally = 1.1 x lO'^^s"^. Deviations from the pure Bose-Hubbard model are about 
2%. Reprinted with permission from 1211 . 



An example how the parameters of the effective Hamil- 
tonian (|20] | vary as a function of the intensity of the driv- 
ing laser f2 is given in figure|6] where the parameters of the 
atom cavity system are chosen to be 324 = .913, N = 1000, 
A = -.913/20, S = 2OOO//VC/13 and a = gis/lO. Fig- 
ure|6]shows the interactions Ub, Uc and Ubc, the tunnelling 
rates Ji,b, Jcc and Jbc and \ fj.c~ l^b\ as a function of fi/gi^. 
For 5 « i? one has \Ubc\ < \Ub\, \Uc\ and Jbb « Jcc ~ 
Jbc- Whenever \fj.c ^ fJ-b\ < \Jbc\, polaritons get con- 
verted into polaritons and vice versa via the tunnelling 
Jbc- With the present choice of a and S, this happens for 
O.I65 < (2 < 1.65. 



4.2. Frozen light: photonic Bose-Hubbard 
models 




10^ 10^ 10-' 



Figure 6 Left: The polariton interactions Ut (dashed line), Uc 
(grey line) and Ubc (solid line) in units of gis as a function 
of n/gis. Right: The tunnelling rates \ Jbb\ (dashed line), \ Jcc\ 
(grey line) and \Jbc\ (solid line) together with |/ic — fib\ (dot- 
ted line) in units of gi3 as a function of H/gi^. The parame- 
ters of the system are (724 = gn, N — 1000, A = —313/20, 
S = 2000\/iV(;i3 and a — gis/lO. Reprinted with permission 
from (381. 



In the last section we have discussed that, as soon as a 
strong coupling regime with a cooperativity factor larger 



than unity, f 3> 1 is achieved, it is possible to engineer po- 
laritonic Bose-Hubbard models in coupled cavity arrays. 
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The prospect that even higher cooperativity factors (of the 
order of 10^ — lO'^ ED) can be realised in some cav- 
ity QED setups, opens up the possibility of engineering 
a Bose-Hubbard model also for the cavity mode photons. 
Such a setting is actually conceptually simpler than the po- 
lari tonic and has already been outlined before. Given that 
the photonic tunnelling term is a consequence of the inter- 
action of neighbouring cavities, to realise a photonic Bose- 
Hubbard model one needs to find a way to efficiently cre- 
ate large photon Kerr nonlinearities using the interaction 
of the atoms with the cavity mode. In this section we re- 
view two proposals to this aim. 

To create a pure photonic Bose-Hubbard model would 
be very interesting as in this way new regimes for light, 
which do not occur naturally, could be engineered. For ex- 
ample, in the photonic Mott insulator exactly one photon 
exists in each cavity, provided the whole structure contains 
on average one photon per cavity. Moreover, the photons 
are localised in the cavity they are in and are not able to 
hop between different cavities. In such a situation pho- 
tons behave as strongly correlated particles that are each 
"frozen" to their lattice site, a system that would corre- 
spond to a crystal formed by light. 

The first proposal for generating strong Kerr photon 
nonlinearities, which we discuss here, is the EIT based set- 
ting analysed in the last section. As shown in ||251 . in the 
regime characterised hy g ^ f2 and ^249 ^ l^-!^!, the 
dark sate polariton po becomes a photon and the nonlin- 
earity, c.f. equation (flST l, is given by 



a 

0.998 r 



U : 



.924.9 g 
\An\ n 



(22) 



The effective decay rate for the photons subject to the non- 
linearity is given by 




c 

0.02 



-0.01 
-0.02 

-0.03 



0.5 


1.0 1.5 


2 




in 10-'^ s 








II 
1 1 










1 1 










0.4 

0.2 



1^-^ ^ . 



(23) 
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Figure 7 Mott insulator state for three photons in three cavi- 
ties: a: photon number ni, b: photon number fluctuation Fi, c: 
differences between full and effective Bose-Hubbard model, Sni 
(solid line) and 5Fi (dashed line) and d: ni and Fi as given by 
a solution of a master equation for the three site Bose-Hubbard 
model (O including damping F. Reprinted with permission from 

m 



Figure 8 The level struc- 
ture for the generation of the 
effective Hamiltonian i26\ . 
ig, Ai) and {A, A2) are 
the Rabi frequencies and de- 
tunings of the cavity-atom 
and laser-atom interactions. 
Reprinted with permission 
from 1261. 




where cavity decay is the main source of loss since 9245 <C 

\An\. 

Taking, for example, the case of toroidal micro-cavities, 
for which the achievable parameters are predicted to be 
524 = 2.5 X lQ^s-\ 74 = 1.6 X 10^s-\ and k = 0.4 x 
10^ ISTI . and setting the parameters such that g/Sl = 
g24g/\Afl\ = 0.1, leads to U/T = 625, which would be 
sufficient to observe a Mott insulator of photons, as shown 
in figure |7] 

A second proposal to realise large Kerr nonlinearities 
is due to Ref. 1261 . There a new configuration, involving a 
three levels lambda structure and only one coupling to the 
cavity mode (see figure [8]l, was shown to produce nonlin- 
earities as large as the EIT based setting discussed before. 
A possible advantage of this approach is that, as the only 
three levels are used, it might be easier to implement it in 
some of the new cavity QED systems, such as quantum 
dots in photonic crystal cavities or Cooper pair boxes in 
stripline cavities. 

The relevant atomic level structure, depicted in figure 
[8] is a system with two metastable states and an excited 



state. The cavity mode couples dispersively only to the — 
2 transition and the levels and 1 are coupled via a far- 
detuned Raman transition. Finally, the detuning associated 
with the lasers and with the cavity mode are assumed to 
be sufficiently different from each other. The origin of the 
nonlinearity can be understood entirely by considering the 
Stark shifts due to the dispersive interactions. 

The full Hamiltonian of the system, in the interaction 
picture, reads 



H = ge-''^^^aS2Q + V2Ae 



-iA2t 



S2+ + h.c. 



Where ^20 := J\ \2k){Qk\ and ^2+ := Ek 



(24) 



As 



shown in figure [8J (g, Ai) and {A, A2) are the Rabi fre- 
quencies and detunings of the cavity-atom and laser-atom 
interactions, respectively. In the dispersive regime, charac- 
terised by < 1 and -^^^ < 1, where furthermore 

VNg, VnA <^ \A2 — Ai\, the processes driven by the 
cavity-atom and laser-atom interactions can be treated in- 
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dependently and the excited state will hardly be populated. 
Hence an effective Hamiltonian for the two metastable states 
and the cavity mode can be found by adiabatically elimi- 
nating the excited state. Dropping out terms proportional 
to the identity, it is given by 



Hi = ——a^aSi 







(5*10 + Sqi), 



(25) 



with e := 2A^/A2. 

In a second interaction picture with respect to Hq — 

^a^'a it reads = ^a^a (5'+_ + 5_+) + f S3, 

where := J^, \+k){-l\, = (5+-)t, and, as 

before, ^3 Ek=i \+k){+k\ - \-k){-kl In the {|±)} 
basis the system can thus be viewed as an ensemble of 
two level atoms driven by a laser with a photon-number- 
dependent Rabi frequency. Considering the dispersive regime 

of this system, <C 1, the atoms prepared in the |— ) 

state will experience a Stark shift proportional to (ata)^, 
which gives rise to a strong Kerr nonlinearity. The effec- 

4 

tive Hamiltonian is be given by Heff — ^^i^ a^aa^aSa. 

Therefore, if all atoms are prepared in the |— ) state, one 
obtains an essentially absorption free Kerr nonlinearity given 
by 



-H^kerr 



2AiO 2Z\] 



-ga^aa^a. 



(26) 



4.3. Measurements 

In this section we discuss ways how information about 
the state of the polaritonic Bose-Hubbard models can be 
obtained in experiments. The possibilities to access the 
system in measurements are somewhat complementary to 
those for cold atoms in optical lattices JH]. Whereas time of 
flight images give access to global quantities such as phase 
coherence in optical lattices, coupled cavity arrays allow 
to access the local particle statistics in measurements. Un- 
like optical lattices they would thus allow to obtain direct 
experimental evidence for the exact integer particle num- 
ber in each lattice site in a Mott Insulator regime. We first 
discuss the measurement scheme for the one-component 
model and then turn to its extension to for the two-component 
model. 



4.3.1. One-component model 

For the model Q, the number of polaritons in one individ- 
ual cavity can be measured via state selective resonance 
fluorescence. To that end, the polaritons are transformed 
into purely atomic excitations by adiabatically switching 
off the driving laser The process is adiabaticif g_B^^^i7 <C 
1/^+1, Im- I 1761 , which means it can be fast enough to pre- 
vent polaritons from hopping between cavities during the 
switching. Hence, in each cavity, the final number of atomic 



excitations in level 2 is equal to the initial number of dark 
state polaritons, while the other polariton species will not 
be transferred to the atomic level 2 (c.f. (fTOl)). In this state 
cavity loss is eliminated and the atomic excitations are 
very long lived. They can then be measured with high ef- 
ficiency and precision using state-selective resonance flu- 
orescence as proposed in ll86l and lISTl . 



4.3.2. Two-component model 



The number statistics for both polariton species for the 
model ( l20t . and c^, in one cavity can again be mea- 
sured using state selective resonance fluorescence. In the 
two-component case the STIRAP can however not be ap- 
plied as in the single component case because the energies 
/ib and are similar and the passage would thus need to 
be extremely slow to be adiabatic. 

For two components, one can do the measurements as 
follows. First the external driving laser i7 is switched off. 
Then the roles of atomic levels 1 and 2 are interchanged in 
each atom via a Raman transition by applying a 7r/2-pulse. 
To this end the transitions 1^3 and 2 ^ 3 are driven 
with two lasers (both have the same Rabi frequency A) in 
two-photon resonance for a time T = 7r(5yi/|vlp {Sa is the 
detuning from atomic level 3). The configuration is shown 
in figure|9h. This pulse results in the mapping | !_,) ^ \2j) 
for all atoms j. 

Next another laser, 0, that drives the transition 1^4 
is switched on, see figure |9j3. Together with the coupling 
(724, this configuration can be described in terms of three 
polaritons, gj, and gl, in an analogous way to pj, p\_ 

and where now the roles of the atomic levels 1 and 
2 and the levels 3 and 4 are interchanged. Hence, if one 
chooses = Q the 7r/2-pulse maps the 6^ onto the dark 
state polaritons of the new configuration, q}^, whereas for 
= — i7 it maps the onto q\. The driving laser is then 
adiabatically switched off, 6* ^ 0, and the corresponding 
STIRAP process maps the q}^ completely onto atomic exci- 
tations of level 1. This process can be fast since the detun- 
ing A is significantly smaller than 5 and hence the energies 
of all polariton species q\, q\_ and q^_ are well separated. 
Another 7r/2-pulse finally maps the excitations of level 1 
onto excitations of level 2, which can be measured by state 
selective resonance fluorescence in the same way as for the 
one-component model. 

The whole sequence of 7r/2-pulse, STIRAP process 
and another 7r/2-pulse can be done much faster than the 
timescale set by the dynamics of the Hamiltonian (|20] | and 
or can be mapped onto atomic excitations in a time 
in which they are not able to move between sites. The pro- 
cedure thus allows to measure the instantaneous local par- 
ticle statistics of each species separately. 
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Figure 9 a: Configuration of the 7r/2-pulse. Two driving lasers 
in two-photon transition with identical Rabi frequencies A couple 
to the atomic transitions 1^3 and 2 <^ 3. b: Configuration 
for the STIRAP process. A driving laser couples to the 1^4 
transition with Rabi frequency 0. The cavity mode couples to 
transitions 2^4 and 1 <-> 3, where the coupling to 1 <-> 3 is 
ineffective and not shown. Reprinted with permission from 1381 . 



5. Models Based on the Jaynes-Cummings 
Hamiltonian 

Many-body dynamics on a lattice that can either be domi- 
nated by the tunnelling between adjacent lattice sites or by 
an on-site interaction can also be generated in an array of 
coupled cavities with Jaynes-Cummings type coupling be- 
tween cavity photons and two level atoms in each cavity. 
As shown in ll22ll and 1231 . the anharmonicity of the energy 
levels of the dressed states of a Jaynes-Cummings Hamil- 
tonian acts as an effective on-site interaction. To see this let 
us recall the Jaynes-Cummings Hamiltonian for one atom, 

H'''' =u;ca^a + LOo\e){e\+gia^g){e\ + a\e){g\), (27) 

where luc , the frequencies of the resonant mode of 

the cavity and of the atomic transition, respectively, g is 
Jaynes-Cummings coupling between the cavity mode and 
the two level system, is the creation operator of a photon 
in the resonant cavity mode, and \g) , |e) are the ground and 
excited states of the two level system. 

Since the energy of a photon ojc and the atomic tran- 
sition energy luq are much greater than the atom photon 
coupling g, the number of excitations is conserved for the 
Hamiltonian ( l27b . Hence it can be diagonalised for each 
manifold with a fixed number of excitations n separately. 
The energy eigenvalues for n excitations read Eq = and 



Et =ncoc + j±]Jng^ + — (28) 

for n> 1, where A = ujq — ujc- 

Let us now consider an array of cavities where the 
atom photon interaction in each cavity is described by equa- 
tion ( [27] ) and assume that we have prepared the system in 
a state with one excitation of energy E^ in each cavity. 
Since the lowest energy for two excitations in one cavity is 
E2 , moving one additional excitation to a cavity requires 
an extra energy of 

E2 = 2^g2 + _ _ ^2.92 + — - y , (29) 



Figure 10 Dependence of the effective repulsion iSj — 2i5j on 
the detuning A. The effective repulsion vanishes for large A 



which plays the role of an effective on-site repulsion. This 
repulsion can be tuned via the detuning A as figure [TOl 
shows. 

The full effective many-body model in a cavity array 
that is based on the effective interaction (|29] l reads, 

^ = E («k«R'+h-c.), (30) 

R <R.R'> 

where R labels the site of a cavity. Each cavity contains 
one atom interacting via the Jaynes-Cummings interaction 
with the cavity mode and photons tunnel between neigh- 
bouring cavities at a rate J. 

Figure[TT]gives evidence for the Mott insulator and su- 
perfluid regimes of Hamiltonian (l30l l. The probabilities to 
find one or two excitations in one cavity are shown for an 
array where initially one cavity had no and all others one 
excitation. For detuning A = Q, double occupations never 
occur, indicating a Mott insulator behaviour. For detunings 
A^ g,on the other hand, occupation numbers larger than 
one occur as to be expected for a superfluid regime. In the 
deep Mott insulator regime, the model (|30] | can also repro- 
duce effective spin models [?], see section ItTI 

Effective many-body physics based on the Hamilto- 
nian (|30] | can not only be observed for a single two level 
system in each cavity, as assumed in (|27] |. but also for 
setups with several two level systems per cavity. Such a 
model can describe photonic crystal micro-cavities doped 
with substitutional donor or acceptor impurities. This ap- 
proach to implementing effective many -body models, which 
can have suitable parameters, has been proposed in |44). 
The phase transitions of a model with several two level 
systems in each cavity have also been studied in P9II47I . 

To allow for the observation of coherent effective many- 
body dynamics, the effective on-site repulsion (|29] | has to 
be much stronger than the rates of spontaneous emission 
and cavity decay. Since the dressed states are superposi- 
tions of a photonic and an atomic component, they are 
vulnerable to both decay mechanisms. Spontaneous emis- 
sion can be suppressed by increasing the detuning A as 
this causes the occupation of the excited atomic level to 
decrease. An increase of A however also causes the re- 
pulsion to decrease as shown in figure [TO] The effective 
many-body dynamics is thus only observable if each cav- 
ity operates in the strong coupling regime, see section|2] 
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Figure 11 Probabilities for finding one (thick solid line) or two 
(thick dashed line) excitations in the middle cavity which was 
initially taken to be empty (in the resonant regime). In this case, 
as there were initially fewer excitations than cavities, oscillations 
occur for the single excitation. The double one is never occupied 
due to the blockade effect-Mott insulator phase (thick dashed 
line). For the detuned case however, hopping of more than one 
excitation is allowed (superfluid regime) which is evident as the 
higher excitation manifolds are being populated now (dashed line 
and thin solid line). Reprinted with permission from 1221 . 



As for the polaritonic Bose-Hubbard model (section 
14.11 ). the transition between Mott insulator and superfluid 
phase can be characterised by the fluctuations of the num- 
ber of excitations in one single cavity. Here, the number of 
excitations in cavity / is measured by the operator J\fi = 
alai + |e;)(e;| and their fluctuations Fi are given by the 
ground state variance of Mi, Fi = {Ni) — {■Mi)'^. Figure 
[T2] shows the number fluctuations Fi, which play the role 
of an order parameter, as a function of the atom cavity de- 
tuning A. Whereas the fluctuations are very small in the 
Mott insulator regime (small detuning), they become sig- 
nificantly larger in the superfluid phase (large detuning). 

Various approaches have shown that this model indeed 
exhibits a quantum phase transition that closely resem- 
bles the Mott insulator to superfluid transition of the Bose- 
Hubbard model Il22ll23]l39ll44ll40ll42l. see also section |6] 
On the other hand deviations of the phase transitions in the 
model ( l30b from the Mott insulator to superfluid transition 
in a Bose-Hubbard model have also been found. 

In this analysis exact diagonalisations of small mod- 
els of very few cavities have been performed BTI . Here, 
insulator states have been found, where the number of ex- 
citations in one site does not fluctuate and the effective 
particles are either atomic (for A/ g < —1 and \ J\ < \A\) 
or polaritonic {A/g < 1 and \ J\ < \A\). Superfluid states 
on the other hand show local excitation number fluctua- 
tions and can either contain polaritonic {A/g < — 1 and 
I J| « \A\) or photonic (A < 0, \J\ > \A\ > g or: A > 
0,\J\ > g) effective particles. 




Figure 12 The order parameter Fi as a function of the detun- 
ing between the hopping photon and the doped two-level sys- 
tem in logarithmic units of the matter-light coupling g. Simu- 
lations include results for 37 sites, with and without dissipation 
due to spontaneous emission and cavity leakage. The graphs are 
obtained by diagonalisation of the Hamiltonian i30\ . 



For the characterisation of the various phases, possible 
signatures in the visibility of interference fringes of pho- 
tons that could be emitted by the cavities have been dis- 
cussed II39II50I . The visibility V of the interference pattern 
is given in terms of the photon number distribution in mo- 
mentum space, S{k) = i Ej^i=i e2'''(^"')/^(a]a;), (L is 
the number of sites in a one dimensional array) by V = 
(max(S') — min(S')) / (max(S') +min(S')). In the super- 
fluid regime, V approaches unity whereas it is very small 
for a Mott insulator 



6. Phase Diagrams 

As already indicated above one of the key questions for 
effective many-body systems in coupled cavity arrays is 
which quantum phase transitions these can show. The phase 
diagram of both, polaritonic Bose-Hubbard models and 
Jaynes-Cummings interaction based models, has been anal- 
ysed in various settings. Approaches based on different 
techniques have been carried out for one- and two-dimensional 
models. 

The phase boundary between a Mott insulator and a 
superfluid phase can be determined in a grand canonical 
approach. To this end, a chemical potential /^i is introduced 
and the grand canonical Hamiltonian reads, 

jjGC = H - fiAf, (31) 

where the operator Af counts the total number of excita- 
tions and is given by 

AA = ^(aJ,aR + |eR)(eR|) (32) 

R 

for Jaynes-Cummings interaction based models (l30l l and 

^f = Y. + |2r> (2rI + |3r> (3r| + 2|4r) (4r|) 

R 

(33) 
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for EIT based models (|9]l, which lead to polaritonic Bose- 
Hubbard models. Here, an excitation of the atomic level 
4 is equivalent to two polariton excitations since an atom 
needs to absorb two cavity photons to make a transition 
from level 1 to level 4, whereas the transition from level 1 
to 2 or 3 only requires one cavity photon. The occupations 
of level 4 thus carry a double weight. 

We note that the chemical potential /i introduced in 
equation dJTT i is conceptually different from the chemical 
potential e in equation ( fTSI l. The approach in section \4~T\ 
considers a closed system without particle exchange with 
the environment and e is an energy shift that affects the 
polaritons. The grand canonical approach in this section, in 
contrast, considers a situation in which particle exchange 
with the surrounding is permitted and is used because of 
its convenience for determining the phase diagrams. 

The total number of excitations is conserved for both 
models, [H, Af] = 0. The boundary between Mott insula- 
tor and superfluid phases is determined by the value of fj, 
for which adding or removing a particle does not require 
energy. 
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Figure 13 Upper panels: Phase diagram for the Hamiltonian 
model ( I30t , with A'^ = 1, 2 atoms inside each cavity at e = cj. 
Lower panels: System compressibility k for the first lobe (i.e., 
p = 1), for different system sizes L, with A'^ = 1 (left) and 
N = 2 (right). Reprinted with permission from 1391 . 



6.1. Phase diagrams for one dimensional 
models 

The phase diagrams for the one dimensional versions of 
both, polaritonic Bose-Hubbard and Jaynes-Cummings in- 
teraction based models, have been obtained by DMRG cal- 
culations (see |[88l for a review on DMRG). These results 
were presented in 1391 . where in addition the existence of 
a polaritonic glassy phase in setups with disorder has been 
predicted. 

The authors consider a linear array of L cavities with 
various numbers of atoms, N, in each cavity. Figure [13] 
shows the phase diagram and compressibility for the model 
(l30b . Here oj denotes the frequency of the cavity-photons 
(u! = oJc), e the energy of the excited atomic levels (e = 
ujo), (3 the atom-photon coupling (/3 = g) and p the number 
of excitations per cavity. The compressibility k is defined 
by K = dp/ dp,. The Mott insulator is characterised by a 
vanishing compressibility whereas k is finite for a super- 
fluid. Figure [14] shows the phase diagram and compress- 
ibility for the model (O. 

For a Bose-Hubbard Hamiltonian, the lobe width at 
t = is given by the on-site repulsion U and the criti- 
cal value t* of the hopping strength at which the various 
lobes shrink into a point is proportional to U. In contrast to 
the Bose-Hubbard model, in both models, ( l30l l and ((9]l, for 
fixed N, the critical values t* are not proportional to the 
lobe width at i = 0. Furthermore, the ratio between the 
upper and the lower slopes of the lobes at small hopping is 
greater than the one predicted for the Bose-Hubbard model 
(2npoiariton for the lower and -2(?7.poiariton + 1) for the up- 
per slope ||89|); this discrepancy disappears on increasing 
the number of atoms inside the cavity. In terms of an ef- 
fective Bose-Hubbard model, this may be interpreted as 
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Figure 14 Upper panels: Phase diagram of model ([9]l. The de- 
tuning parameters are set to zero and g/Q — 1. Lower panels: 
System compressibility k for the first lobe, for different system 
sizes L, with A'^ = 1 (left) and N — 2 (right). Reprinted with 
permission from 1391 . 



a correlated hopping of the polaritons, where the hopping 
depends on the occupation of the cavity. For larger occu- 
pations, the polaritons have a stronger tendency to tunnel 
between cavities and consequently the superfluid regime is 
already reached for lower values of the hopping term pref- 
actor t. Deviations of the critical hopping t* from the crit- 
ical hopping in a Bose-Hubbard model vanish for > 10 
atoms per cavity 1391 . 
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For the detection of the Mott insulator-superfluid quan- 
tum phase transition, the polariton and photon fluctuations 
in a block of linear dimension AI (in units of the lattice 
constant of the array) have been studied P9ll . Whereas po- 
lariton fluctuations are independent of the block size M in 
the Mott insulator regime, they grow logarithmically with 
M in the superfluid phase. The prefactor of the logarithm 
is thereby related to the compressibility of the system. For 
the photon fluctuations, the scaling of the leading terms is 
always linear in M and non-critical while the critical be- 
haviour is encoded in the subleading scaling P9l . 



6.2. Phase diagrams for two dimensional 
models 

The phase diagrams for the two dimensional version of the 
Jaynes-Cummings based model ( l30l l have been calculated 
by various approaches, among those a mean field approach 
and a variational cluster approach [i40l . 



6.2.1. Mean-field phase diagram for the 2d 
Jaynes-Cummings based model 

A mean-field phase diagram for the two dimensional ver- 
sion of the model ( l30l l on a triangular lattice has been 
obtained in 1231 . Introducing the superfluid order param- 
eter <F = (or) and using the decoupling approximation 

aROR' = (aR)aR' + aR(aR') + (ak)(»R.')' obtains 
the mean-field Hamiltonian, 



H 



MF 



R 



zJW (aL + aR ) + zJltf- 



M (aRflR + |eR)(eR|) , 



(34) 



where z counts the number of nearest neighbours and is 
here set to z = 3 since a triangular lattice is assumed 1231 . 
Here again, a chemical potential jj, has been introduced in 
(|34] | and the phase boundary between superfluid and Mott 
insulator is given by the value of /i, for which adding or 
removing one excitation to or from the system does not 
cost energy. 

The mean-field phase diagram for the two dimensional 
model ( |34] | on a triangular lattice is shown in figure [15] 
Here, 'F = corresponds to stable Mott lobes, where the 
number of excitations in each cavity increases with fi. The 
regions on the right correspond to ^ 0, where the sys- 
tem will be found in a superfluid phase, i.e. within the 
mean field approximation, the ground state corresponds to 
a coherent state of photons in each cavity. The number of 
photons in each Mott lobe can be confirmed by consider- 
ing the average number of excitations per site in the grand 
canonical ensemble p = dEg/dfi, where Eg is the ground 
state energy. Figure [16] shows p as a function of the hop- 
ping and chemical potential. The plateaus indicate regions 
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Figure 15 Slices showing the superfluid order parameter, •f', as 
a function of the photon hopping frequency, here denoted by n 
(k = J), and the chemical potential, /i. For detuning, A — (a), 
A — 2/3 (f3 denotes the Jaynes-Cummings coupling, (3 = g) (b) 
and A = —2j3 (c) (Note that the figure defines A = ujc — uJo 
and denotes the cavity resonance by to, contrary to the text). Mott 
insulator lobes are indicated by the regions of9 = Q, the low- 
est of which are shown. Dominating the left-hand edge (where 
photonic repulsion dominates over hopping) is the Mott insu- 
lator phase (denoted Ml), and the superfluid phase is found on 
the right-hand edge (denoted SF). The white contour in a corre- 
sponds to the region where <!> becomes non-zero, delineating the 
quantum phase transition. Reprinted with permission from 1231 . 




Figure 16 Plateaus with constant density, p, indicating regions 
with a definite number state excitations, as a function of the 
chemical potential p and photon hopping frequency k for detun- 
ing A = Q. Regions with varying p have coherent states as the 
ground state configuration. Reprinted with permission from 1231 . 
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t/g 

Figure 17 Ground-state phase diagram for the Hamiltonian OOt : 
(a) ID, (b) 2D and (c) different coordination numbers z. Lines are 
guides to the eye. Reprinted with permission from |40|. 



of constant density and incompressibility which are char- 
acteristic for a Mott insulator phase. 

A mean-field phase diagram for settings with more than 
one two level system in each cavity has been calculated in 
ll44l . where 8 two level systems par cavity have been con- 
sidered. Such a model describes photonic crystal micro- 
cavities doped with substitutional donor or acceptor impu- 
rities. 

6.2.2. Phase diagrams for 2d models from a variational 
cluster approach 

The phase diagram of model dSOl l for the one and two di- 
mensional case and the stability of its Mott insulator states 
with respect to non-zero temperatures have been analysed 
in ||40| . The authors employ a variational cluster approach 
that is expected to become exact with increasing cluster 
size Lc. FigurefTTlshows the phase diagram in one dimen- 
sion, two dimensions and for different coordination num- 
bers z. fi is the chemical potential, g the Jaynes-Cummings 
coupling and t the tunnelling rate of photons between neigh- 
bouring cavities. Figure[T8]shows the unit filling Mott lobes 
for both, the one and two dimensional model ( l30l l at finite 
temperatures T. 
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Figure 18 Mott-like regions for different T for the Hamiltonian 
i30l . Mott lobes shrink with increasing temperature T. Reprinted 
with permission from go). 



Phase diagrams for the model ( [30] l in two dimensions 
have also been obtained by the stochastic series expansion 
quantum Monte Carlo method ll43l . where it has been pre- 
dicted that the universality class of coupled cavity arrays 
may be different from the Bose-Hubbard model. 



6.3. Bose-glass phases 

A glassy phase for polaritons in a one dimensional cavity 
array with disorder has been predicted in ll39l . This phase 
is characterised by a finite compressibility, gapless exci- 
tation spectrum and zero superfluid density and has been 
predicted to exist in a Bose-Hubbard model for relative 
fluctuations of the interaction, SUcff/ (C4ff), uniformly dis- 
tributed in an interval of length 2e = 0.5 in a system of 
L = 200 sites at filling p ^ 1.01 for 0.078 < t^ff/ (f/off) < 
0.133. Figure [19] shows relative fluctuations of the interac- 
tion SUes/ (C/off) as a function of atom number fluctua- 
tions over an interval SN in the cavities. For the model 
( l30l l. a glassy phase should be visible for , e.g. (N) = 100 
particles per cavity and SN ~ 19 (such to have e = 0.25) 
in the range 7.51 x lO'^ < t//3 < 1.401 x lO'^ . 

The observability of quantum phase transitions from 
a superfluid to a Bose-glass to a Mott insulator that is 
driven by an decrease of the photon tunnelling rate J has 
also been predicted for a two dimensional model in B4l . 
Here, photonic crystal micro-cavities doped with substi- 
tutional donor or acceptor impurities have been consid- 
ered, where fluctuations in the number of impurities (AN), 
the impurity-photon coupling (Agk) and the detunings be- 
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Figure 19 Relative standard deviation of the effective on-site 
interaction strength UcS averaged over L — 10* cavities, as a 
function of the atom number fluctuations SN; from left to right 
(iV) = 1, 10, 100, 1000. Dashed line indicates an effective 
variation of the interaction strength equal to the standard devia- 
tion of a random interaction uniformly distributed in the interval 
Uesii)/ (f/eff) G [— e, +e] with e = 0.25. In the inset an exam- 
ple of the variation on the on-site effective interaction is shown; 
{N) = 100, 5N = 20. Reprinted with permission from OH. 




< 
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0.2 I ■ 0.2 

AN/V<N> 0.3 0,3 Ag^/g 

Figure 20 The isovalue surface that separates the existence and 
nonexistence of an equilibrium polaritonic quantum phase transi- 
tion from the Bose glass to Mott insulator states. Reprinted with 
permission from 1441 . 



tween impurity and photon resonances (Z\ijjph) are expected 
in experimental realisations. The information whether an 
equilibrium polaritonic quantum phase transition from a 
Bose glass to Mott insulator states exists is obtained by 
plotting the isosurface in which the critical polariton tun- 
nelling rate at the Bose glass Mott insulator phase bound- 
ary is set to be 10 times larger than the polariton loss rate, 
assuming a cavity Q-factor 10^, an impurity number N = 
3 and an impurity photon detuning A = 12g. Inside the 
isovalue surface, quantum phase transitions from a super- 
fluid to a Bose glass to Mott insulator states can be ob- 
served by continuously decreasing the tunnelling rate J. 



Outside the isovalue surface, only superfluid and Bose glass 
phases exist. 

7. Effective Spin Models in Coupled Cavity 
Arrays 

In the previous sections we have seen that under appro- 
priate conditions interesting bosonic many-body Hamilto- 
nians can be created and probed in coupled cavity arrays. 
There the atoms were used for detection and most impor- 
tantly to generate an interaction between photons in the 
same cavity, either actively, as in the polaritonic case, or in 
an indirect manner, as for the creation of large Kerr non- 
linearities. In this section we focus on a complementary 
regime, in which the tunnelling of photons is used to me- 
diate interactions for atoms in neighbouring cavities. We 
review proposals to realise effective spin lattice Hamilto- 
nians with individual atoms in micro-cavities that are cou- 
pled to each other ll27ll29ll3T1 . 

Interacting spin systems have a key role in quantum 
statistical mechanics, condensed-matter physics, and more 
recently in quantum information science. On the one hand 
they have a simple microscopic description, which allows 
the use of techniques from statistical mechanics in their 
analysis. On the other, they have a very rich structure 1901 
and can hence be employed to model and learn about phe- 
nomena of real condensed-matter systems. In quantum in- 
formation science, spin chains can be harnessed to prop- 
agate and manipulate quantum information, as well as to 
implement quantum computation. We first discuss how spin 
models that conserve the total number of excitations can 
be generated as limiting cases of Bose-Hubbard type mod- 
els. Then we review an approach where the full anisotropic 
Heisenberg Hamiltonian can be engineered for atoms that 
interact via the exchange of virtual photons. Finally we 
mention schemes that investigate magnetic properties for 
photons of different polarisations. 

7.1. Spin models that conserve the total number 
of excitations 

A certain class of effective spin Hamiltonians can be gen- 
erated as a limiting case of the effective many -body Hamil- 
tonians discussed in sections 2] and |5] Whenever the engi- 
neered on-site interactions are repulsive and strong enough 
to suppress the occupation of states with more than one 
excitation almost completely, the states with zero and one 
excitation can be identified with | f) and | [) respectively 
II22I and the dynamics can effectively be described by the 
Hamiltonian 

H-T.B-^+ E (35) 

j <j,j'> 

where the second sum runs over nearest neighbours. ctJ^^^^ 
are the Pauli matrices at site j, J a coupling coefficient and 
B a magnetic field. 
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The Hamiltonian i35[ conserves the number of spins 
in the state | f) and can thus be solved efficiently. Much 
richer spin Hamiltonians can be generated when putting 
atoms of a A level structure into the cavities. 



7.2. The anisotropic Heisenberg model 



A particular interesting spin Hamiltonian is the anisotropic 
Heisenberg or XYZ model, given by 

<j-r> 

This model contains a rich phase diagram and can be used 
to study quantum magnetism as a limiting case of the fermic 
Hubbard model, which is believed to contains the main 
features of high Tc superconductors ||90| . From a quan- 
tum information perspective, such a model can be used to 
create cluster states, highly entangled multipartite quan- 
tum states which are universal for quantum computation 
by single-site measurements ||9T| . 

The Hamiltonian ( [36l ) can be generated in an array of 
coupled cavities where each cavity contains one atom of a 
A level structure, c.f. figure|2T|and[22] The two spin polar- 
isations I t) and I I) are represented by the two longlived 
atomic levels of the A level structures. Together with exter- 
nal lasers, the cavity mode that couples to the atom inside 
each cavity can induce Raman transitions between these 
two longlived levels. Due to a large enough detuning be- 
tween laser and cavity mode, these transitions can only 
create virtual photons in the cavity mode which mediate 
an interaction with another atom in a neighbouring cav- 
ity. With appropriately chosen detunings, both the excited 
atomic levels and photon states have vanishing occupation 
and can be eliminated from the description. As a result, the 
dynamics is confined to only two states per atom, the long- 
lived levels, and can be described by a spin 1/2 Hamilto- 
nian. The small occupation of photon states and excited 
atomic levels also strongly suppresses spontaneous emis- 
sion and cavity decay. As discussed later, this proposal can 
be realised as soon as the strong coupling regime of cavity 
QED is achieved. 

We note that similar ideas have been recently proposed 
to realise spin systems of larger spins |f3Tl and with inter- 
actions that depend on the direction of the respective lattice 
edge 1291. 

In the next subsections we first discuss how effective 
(T^cr^, a^cry and a^cr^ interactions together with an ef- 
fective magnetic field Baz can be engineered and then 
explain how these can be combined to generate the full 
anisotropic Heisenberg model ( |36] |. 




a 



Figure 21 Level structure, driving lasers and relevant cou- 
plings to the cavity mode to generate effective a^a^- and cr^a^- 
couplings for one atom. The cavity mode couples with strengths 
Qa and Qb to transitions \a) ^ |e) and \b) <-» |e) respectively. 
One laser with frequency tUa couples to transition \a) ^ |e) 
with Rabi frequency Oa and another laser with frequency ujb to 
1 6) ^ je) with Ob- The dominant 2-photon processes are indi- 
cated in faint grey arrows. Reprinted with permission from II27I . 



7.2.1. XX and YY interactions 

XX and YY interactions can be generated in an array of 
coupled cavities with one 3-level atom in each cavity (see 
figure I2TI1. The two long lived levels, \a) and \b), represent 
the two spin states. The cavity mode couples to the tran- 
sitions I a) ^ |e) and \b) ^ |e), where |e) is the excited 
state of the atom. Furthermore, two driving lasers couple 
to the transitions \a) ^ \e) respectively \b) ^ |e). For the 
sake of simplicity we consider here a one-dimensional ar- 
ray, but the generalisation to higher dimensions is straight 
forward. The Hamiltonian of the atoms reads 

N 

where the index j counts the cavities, uj^ is the energy of 
the excited level and coab the energy of level |6). The en- 
ergy of level I a) is set to zero. The Hamiltonian that de- 
scribes the photons in the cavity modes is given by equa- 
tion. For simplicity of notation we set Jc — 2u!ca and 
for convenience we assume periodic boundary conditions, 
where He can be diagonalised via the Fourier transform 

1 .fc, , N ^ , ^ N 

VN ^ N ' 2 - - 2 

i=l 

for N odd to give He = Sfc^fe'^I'^fc '^^^^ ^fe ^ + 
2Jc cos(fc). Finally the interaction between the atoms and 
the photons as well as the driving by the lasers are de- 
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scribed by 

N 

Hac = J2 



'■*+.g„a, )|e,)(a,| (38) 



+ (^e-'-^'+g,a,]\e,){b,\+h.c. 



Here ga and gt, are the couplings of the respective transi- 
tions to the cavity mode, Ha is the Rabi frequency of one 
laser with frequency uJa and f2i, the Rabi frequency of a 
second laser with frequency iUb ll92l . The complete Hamil- 
tonian is then given hy H = Ha + He + Hac ■ 

By switching to an interaction picture with respect to 
Ha = Ha + He - Si \bj){bj\, where Si = ujab - 

{u!a — cJb)/2, the excited atom levels \ej) and the photons 
can be adiabatically eliminated 1931 . To obtain effective 
spin-spin interactions, terms up to 2nd order need to be 
considered while dropping fast oscillating terms. For this 
approach the detunings Aa = uJe — oja, Ai, = uj^ — ij->b — 

{uJab - Si), S^ = UJe - LOk and si = LOe- i^k - {Ulab " 

Si) have to be large compared to the couplings Qa, ^b, ga 
and gt- Furthermore, the parameters must be such that the 
dominant Raman transitions between levels a and b are 
those that involve one laser photon and one cavity photon 
each (c.f. figure I2TI1. To avoid excitations of real photons 



via these transitions one requires | Aa — (5; 

figgb ^bg. 
2A. 



\Au 



s':.\ » 



2Au 



Hence whenever the atom emits or absorbs a virtual 
photon into or from the cavity mode, it does a transition 
from level \a) to \b) or vice versa. If one atom emits a vir- 
tual photon in such a process that is absorbed by a neigh- 
bouring atom, which then also does a transition between 
I a) to \b), an effective spin-spin interaction has happened. 
Dropping irrelevant constants, the resulting effective Hamil- 
tonian reads 

N 

H,y = J2 + iJi<J+a-^i+J2(T~(T~+i+h.c.), (39) 



9a 



9b 



- t^ab 



Figure 22 Level structure, driving lasers and relevant couplings 
to the cavity mode to generate effective cr^cr'-couplings for one 
atom. The cavity mode couples with strengths Qa and Qb to tran- 
sitions I a) ^ |e) and \ h) ^ \e) respectively. Two lasers with fre- 
quencies u} and V couple with Rabi frequencies ila respectively 
Aa to transition \a) <-> \e) and respectively /!(, to |6) <-> |e). 
The dominant 2-photon processes are indicated in faint grey ar- 
rows. Reprinted with permission from 1271 . 



For f2a ~ zk{Aaga/Abgb)f2b with Ha and fit, real, 
the interaction is either purely c^cr^ (+) or purely cr^cr^ 
(— ) and the Hamiltonian ( |40] i becomes the Ising model in 
a transverse field, whereas the isotropic XY model {Jx = 
Jy, i.e. J2 = 0) is obtained for either fia 01 f^i, ^ 0. 
The effective magnetic field B in turn can, independently 
of Jx and Jy, be tuned to assume any value between \B\ ^ 
\Jx\, \ Jy\ and \B\ <C \ Jx\, \ Jy\ by varying Si. Hence the 
system can be driven through quantum phase transition by 
varying the Rabi frequencies and detunings of the driving 
lasers. Now we proceed to discuss effective ZZ interac- 
tions. 



where ct| = \bj){bj \ — \aj){aj\, = \bj){aj \ and B, Ji 
and J2 are functions of the various Rabi frequencies and 
detuningfl If = J2, this Hamiltonian reduces to the 
XY model. 



7.2.2. ZZ interactions 



N 



H^y = J2 M + •^-''I'^l+i + ^i^'^'^+i ' (40) 
with Jx = (Ji + J2)/2 and Jy = (Ji - J2)/2. 



^ To second order they are given by B = ^ 



IfiaT 2 2 

2 2 io_|2„2 



and 



1 rifi.p 

J2 = 1^^^^, where 7... = ^Li:i^-7i = 
J. V 1 and ^, = i y siEfiii 



Effective a^a^ interactions can be generated with the same 
atomic level configuration as XX and YY terms, but now 
only one laser with frequency uj mediates atom-atom cou- 
pling via virtual photons. A second laser with frequency 1/ 
is used to tune the effective magnetic field via a Stark shift. 
The atoms together with their couplings to cavity mode 
and lasers are shown in figure |22] Again, we consider the 
one-dimensional case as an example. The generalisation 
to higher dimensions is straightforward. The Hamiltonians 
Ha of the atoms and He of the cavity modes thus have the 
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same form as above, whereas Hac now reads 

N r 



3 = 1 



2 



e-"* + ^e-*+g,a,. )|e,)(a,| 



(41) 



^e--* + ^e-*+5fca,)|e,>(6,|+h.c. 



Here, J?^ and f^i, are the Rabi frequencies of the driving 
laser with frequency oj on transitions \a) \e) and \b) 
|e), whereas Aa and Ai, are the Rabi frequencies of the 
driving laser with frequency on transitions \a) |e) 
and |6) ^ |e). 

In an interaction picture with respect to Hq = Ha + 
He, the excited atom levels \ej) and the photons can be 
adiabatically eliminated ll93l . Again, the detunings Aa = 

LOe - LO, Ab = UJe - UJ - LUab, Aa = UJe - V, Ab = 

uje-v-u!ab,&a = t^e - ^fc and (5^ = ti^g - Ljfc - cJq6 have to 
be large compared to the couplings i?^, J?;,, /1q, /It, Pa and 
gb, whereas now Raman transitions between levels a and 
h should be suppressed. Hence parameters must be such 
that the dominant 2-photon processes are those that involve 
one laser photon and one cavity photon each but where the 
atom does no transition between levels a and h (c.f. fig- 
ure |22l). To avoid excitations of real photons in these pro- 



nag a 

2Aa 



ilbgb 
2Ah 



cesses, one requires \ Aa — <5^| , 

Whenever two atoms exchange a virtual photon in this 
scheme, both experience a Stark shift that depends on the 
state of the partner atom. This conditional Stark shift plays 
the role of an effective tr^cr^ -interaction. Dropping irrele- 
vant constants, the effective Hamiltonian reads 



H„ — 



N 



(42) 



where B and Jz are determined by the detunings and Rabi 
frequencies of the employed laserfl Here again, the in- 
teraction J, and the field B can be tuned independently, 
either by varying Qa and Qb for Jz or by varying Aa and 
Ab for B. In particular, \Aa\'^ and \Ab\'^ can for all values of 
Qa and Qb be chosen such that either Jz ^ B or Jz ^ B. 

Both, XX and Y Y interactions and ZZ interactions can 
be combined to engineer the full Hamiltonian ( l36b . This 
can be done via the Suzuki-Trotter formula. 
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Figure 23 The occupation probability p(ai) of state |ai) (solid 
line) and the probability p{li) of spin 1 to point down (dashed 
line) for the parameters tOe — lO^'GHz, cuab = 30GHz, Aa ~ 
30GHz, At = 60GHz, ujc ^ - At + 2GHz, Aa = 15GHz, 
na = (ib = 2GHz, Aa ^ Ab = 0.71GHz, ga = gb = IGHz, 
Jc = 0.2GHz and 5i = -0.0165GHz (plot a) respectively 5i = 
—0. 0168GHz (plot b). Both, the occupation of the excited atomic 
states (|ej) (ej I) and the photon number (a^a) are always smaller 
than 0.03. Reprinted with permission from |27| . 



7.2.3. The complete effective model 

Making use of the Suzuki-Trotter formula, the two Hamil- 
tonians (|40] | and ( l42b can now be combined to one effec- 
tive Hamiltonian. To this end, the lasers that generate the 
Hamiltonian ( |40| | are turned on for a short time interval 
dt (I jiJxyl I ■ dt <^ V) followed by another time interval dt 
(I |i?zz| I • <C 1) with the lasers that generate the Hamil- 
tonian (|42]) turned on. This sequence is repeated until the 
total time range to be simulated is covered. 

The effective Hamiltonian simulated by this procedure 
is i/spin = Hjiy + iJzz, which is precisely the Heisen- 
berg anisotropic model of equation. ( l36b with Btot ^ B + 
B. The time interval dt should thereby be chosen such 
that Q^^,g-^ < dti,dt2 < J~^ , J~^ , J~^ , B'^ and 
B^^, so that the Trotter sequence concatenates the effec- 
tive Hamiltonians Hxy and Hzz- The procedure can be 
generalised to higher order Trotter formulae or by turning 
on the sets of lasers for time intervals of different length. 

The validity of all above approximations is shown in 
figure |23] where numerical simulations of the dynamics 
generated by the full Hamiltonian H are compared it to 
the dynamics generated by the effective model ( l40l i. The 
present example considers two atoms in two cavities, ini- 
tially in the state -^(lai) + l^i)) <8i |a2), and calculates the 
occupation probability p(ai) of the state |ai) which cor- 
responds to the probability of spin 1 to point down, p{ii)- 
Figure l23h shows p(ai) and p{ii) for an effective Hamil- 
tonian (ESJ with Bft = 0.135MHz, = 0.065MHz, 



Jy = 0.007MHz and Jz = 0.004MHz and hence |Btot| > 
\Jx\, whereas figure l23b shows p(ai) andp(|i) for an ef- 
fective Hamiltonian (|36] | with the same J^, Jy and Jz but 
Btot = -0.025MHz and hence |Bto,| < \ Jx\ mi- 
Discrepancies between numerical results for the full 
and the effective model are due to higher order terms for 
the parameters B, B, J^, Jy and Jz, which lead to relative 
corrections of up to 10% in the considered cases. Despite 
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this lack of accuracy of the second order approximations, 
the effective model is indeed a spin- 1/2 Hamiltonian as 
occupations of excited atomic and photon states are negli- 
gible. 



7.2.4. Experimental requirements 

For an experimental implementation, the parameters of the 
effective Hamiltonian, J^, Jy, Jz, B and B have to be 
much larger than rates for decay mechanisms via the pho- 
tons or the excited states \ej). With the definitions Q ~ 
max(]7a, g = max(ga,gb), ^ = m\\\{Aa,Ab), the 
occupation of the excited levels \ej) and the photon num- 
ber Up can be estimated to be (|ej)(ej |) « |i7/2Z\p and 
Up « |(i7g/2Z\)7ip, whereas the couplings J^, Jy and Jz 
are approximately \ {Qg /2A)\'^^2- 

Spontaneous emission from the levels \ej) at a rate Fe 
and cavity decay of photons at a rate Fc thus lead to ef- 
fective decay rates 

A = \^2/2A\''Fe and F2 - \{[2g/2A)^i\'Fc. 

The effective dynamics of Hamiltonian ([36b can thus be 
observed if 

A « \{Qg/2A)\''-t2 and F2 « \{[2g l2A)\''-i2, 

which implies Fe <C Jcg^l^^ and Fc <C Jc {Jc < 
5/2), where, 5 = | (cj^ + tJ6)/2 - tJcl for the XX and YY 
interactions and 5 =\ijj — loc] for the ZZ interactions and 
the approximations |7i| w 5^^ and I72I ~ JcS~^ have 
been used. Since photons should be more likely to tunnel 
to the next cavity than decay into free space, Fc <C Jc 
should hold in most cases. For Fe ^ Jcg^ /5^, to hold, 
cavities with a high ratio g / Fe are favourable. Since 6 > 
2Jc, the two requirements together imply that the cavities 
should have a high cooperativity factor ^. 



7.3. Magnetic dynamics of polarised light 

In contrast to identifying the internal states of the atoms 
with spin polarisations, magnetic dynamics can also be 
studied by considering cavities which support two distinct 
modes of light with different polarisations, 1 and |. In ll28ll 
the dynamics of the photon polarisations in such a scheme 
was analysed. The authors consider cavities that are doped 
with atoms of a level structure, where two distinct inter- 
nal transitions only couple to photons of different polari- 
sations. The interactions of the photons with the atoms are 
chosen such that double occupancies are suppressed by the 
blockade effects discussed in section |5] and the dynamics 
of the photons with creation operators aj and aj for the 
two polarisations obeys the effective Hamiltonian 1281 

H= J{S^S=^,+Spy) + J'S]S^,, (43) 

<jj'> 




10 20 30 40 50 60 10 20 30 40 50 60 



tinlO^'^s tinlO^'^s 

Figure 24 a: The von Neumann entropy An of the reduced den- 
sity matrix of 1 effective spin in multiples of In 2 and b: the pu- 
rity of the reduced state of the effective spin model for 3 cavities 
where A = 0.021MHz. The plots assume that no spontaneous 
emission took place. For a spontaneous emission rate of O.lMHz 
{g = IGHz), the probability for a decay event in the total time 
range is 1.5%. Hence, cluster state generation fails with proba- 
bility 0.005 X n forn <C 1/0.005 = 200 cavities, irrespectively 
of the lattice dimension. Reprinted with permission from 1271 . 



whereS"! = ^{a^^a^^-al^aji), = ^Ht'^J'i+sVjT)' 
Sj = 5i(a]-fajj, — a^j^Qj-^) and the indices j and j' label 
the cavities. 

7.4. Cluster state generation 

Effective spin models in arrays of coupled cavities can be 
used for the generation of cluster states ||9T1 , which form, 
together with the local addressability, a platform for one- 
way quantum computation [ 94 , 95 1 . The experimental veri- 
fication that the desired entanglement has indeed been cre- 
ated is thereby a non-trivial task in its own right |!96l. 

One way to generate cluster states is via the effective 
Hamiltonian (|42] | 1271 . A similar rout has been presented 
in l30l . To this end, all atoms are initialised in the states 
{\aj) + \bj)) I\f2, which can be done via a STIRAP process 
l76l . and then evolved under the Hamiltonian ( |42] | for t = 
tt/AJz. 

Figure [24] shows the von Neumann entropy of the re- 
duced density matrix of one effective spin An and the 
purity of the reduced density matrix of the effective spin 
chain for a full three cavity model. Since An ~ log2'^ 
for t w 50/is while the state of the effective spin model 
remains highly pure (A = t^ip^] > 0.995) the degree of 
entanglement will be very close to maximal. Thus the lev- 
els \aj) and \bj) have been driven into a state which is, up 
to local unitary rotations, very close to a three-qubit cluster 
state. 

A second way to generate cluster states is to imple- 
ment controlled phase gates between the qubits in differ- 
ent cavities l52l . First all qubits are initialised in the state 
1+) = (|l) + |0))/V2,where|l)and|0)referto states with 
one or zero excitations in a cavity. The phase gates are then 
implemented by isolating chains of three qubits via a Stark 
shift that turns surrounding qubits off resonance. The iso- 
lated three-qubit chains evolve for a time t ^ tt/ (2\/2J) 
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Figure 25 A 2d cavity array. Con- 
secutive application of gates A, B, 
C and D, depicted in parts (a), 
(b), (c) and (d), each isolate chains 
of 3 qubits, realising controlled- 
phases and SWAPs (swapping quan- 
tum/classical gates) between the 
qubits at either end of the chain 
(the dashed lines indicate where the 
controUed-phase gates have been ap- 
plied). Reprinted with permission 
from (52|. 



(J is the photon tunnelling rate) and the central qubit is 
subsequently measured. Local unitaries on the remaining 
two qubits that depend on the outcome of this measure- 
ment complete the controlled phase gate. Figure |25] shows 
the entire sequence of gates. 

Further ways to create entangled states include mea- 
surement based schemes. One version uses two atoms in 
two cavities that interact via the exchange of photons ||55l . 
Both atoms can emit photons of two polarisations (left cir- 
cular and right-circular), depending on which excited state 
they are in. Hence the two atoms are prepared such that 
they each emit a single photon, one with left-circular and 
one with right-circular polarisation. After a time t ~ n/iJ 
(J being the photon tunnelling rate), when photons have 
tunnelled between the cavities with probability 1/2, the 
two atoms are employed to perform a non-demolition mea- 
surement to determine the number of photons in each cav- 
ity. Provided each cavity contains one photon the state is 
collapsed to a two photon polarisation entangled state thus 
providing a heralded source of entangled photons. 

Apart from quantum information motivated applica- 
tions coupled cavity arrays with local addressability can 
also be used to study novel physical setups. Among those 
are interface between regions in different condensed mat- 
ter phases, see section|9l 



8. Promising Platforms for Realisations 

In this section, we present some cavity QED setups which 
have favourable features for the realisation of the effec- 
tive many-body models described in this review. These 
platforms are thus promising candidates for experimental 
demonstrations of the theory presented here. The selection 
we present here obviously reflects our personal view. It is 
of course possible that devices which we do not explicitly 
mention here may in the end provide the first realisation of 
polaritonic many-body systems as discussed in previous 
sections. 

A major part of all cavity QED experiments have been 
performed with Fabry-Perot cavities formed by rather mas- 
sive mirrors. With such miiTors, coupling two cavities such 
that photons can tunnel between them seems quite diffi- 
cult. On the other hand several new cavity QED structures 



with small volume micro-cavities have emerged in recent 
years. These are routinely produced in large arrays and by 
construction can allow for efficient photon tunnelling be- 
tween neighbouring cavities. In addition, the strong cou- 
pling regime with cooperativity factors much larger than 
unity has been achieved in several of these devices. These 
two features together make such micro-cavity arrays very 
promising candidates for the generation of the effective 
many-body models discussed in this review. In the follow- 
ing we discuss some platforms in a bit more detail. 



8.1. Photonic band gap cavities 

One promising device consists of quantum dots coupled to 
photonic band gap defect nanocavities. Photonic crystals 
are structures with periodic dielectric properties which af- 
fect the motion of photons, in a similar manner as the peri- 
odicity of a semiconductor affects the motion of electrons 
1971 . Such structures can exhibit band gaps in frequency 
space in which no photon of particular frequencies can 
propagate in the material (for an introductory review see 

Ell). 

A nanocavity can be created in a photonic-band-gap 
material by producing a localised defect in the structure 
of the crystal, in such a way that light of a particular fre- 
quency cannot propagate outside the defect area. Large ar- 
rays of such nanocavities have been produced |f99l. Photon 
hopping between neighbouring cavities has been observed 
in the microwave and optical domains 1 100 101 1, where its 
description in the form as given in equation ^ was found 
to be accurate. Quantum dots, semiconductors whose exci- 
tations are confined in all three directions, behave in many 
respects like atoms and can be effectively addressed by 
lasers 11021 . Quantum dots can be created inside photonic 
crystal nanocavities II103II and made to interact with the 
cavity mode to form a standard cavity QED system. The 
strong coupling regime with cooperativity factors up to 
^ w 10 has already been achieved in such systems 11611 . 
Cooperativity factors of ^ ~ 10^ are predicted to be re- 
alisable ISll . Due to the extremely small volume of the 
nanocavity, the Jaynes-Cummings coupling coefficient can 
be extremely large 16011 104L while spontaneous emission 
rates are very low. Figure |26] shows a quantum dot in a 
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Figure 26 A quantum dot 
in a photonic band gap cav- 
ity, a, AFM topography of a 
photonic crystal nanocavity 
aligned to a hill of material 
on the surface arising from 
a quantum dot buried 63 nm 
below, b, Electric field in- 
tensity of the photonic crys- 
tal cavity mode showing that 
the location of the buried 
QD, indicated by the teal 
dot, overlaps the field max- 
imum. Reprinted with per- 
mission from 1611. 



photonic band gap cavity. For the applications discussed 
in this review current cavity decay rates of photonic band 
gap cavities remain a limiting factor This is mainly due to 
the fact that most structures employ two dimensional pho- 
tonic crystals as placing a quantum dot into a three dimen- 
sional structure is difficult. The difficulties to control the 
properties of individual quantum dots and to place them at 
desired locations could even pose a more significant prob- 
lem. Current experiments pick a quantum dot from an en- 
semble and then fabricate a photonic crystal tailored to this 
dot around it. Such an approach can obviously not be used 
for producing a whole array of sufficient quality. 

One might think about trapping atoms in photonic band 
gap cavities instead of letting them interact with quantum 
dots. Since the frequencies of currently available photonic 
band gaps do not fit to atomic transition frequencies, such 
a strategy would require the development of new photonic 
band gap materials. 



8.2. Silicon structures: micro discs and micro 
toroids 




Figure 27 SEM picture of a 
toroidal microresonator with 
an intrinsic cavity Q of 10*. 
Reprinted with permission 
from (m. 



atoms are placed close to the cavity surface and interact 
with the evanescent field. This has been seen in an exper- 
iment where atoms fell by the cavity surface 1581 and a 
photon blockade effect generated by one atom has been 
demonstrated 1108 1. Theoretical predictions [81 J show that 
very large cooperativity factors of ^ ~ 10^ are expected to 
be achieved in such systems. 

For realising effective Hamiltonians such as (|7]i, dSOl l 
or (|36] |, some aspects still remain challenging for micro 
disc and micro toroidal cavities. One aspect is that atoms 
need to be trapped for a sufficiently long time close enough 
to the cavity surface to maintain a sufficient strong cou- 
pling regime. A second aspect is that all cavities of an ar- 
ray need to be tuned into resonance with each other with 
enough accuracy. This can be achieved by varying their 
temperature and thus shrinking or expanding them, but is 
still a technical challenge. Finally there can be modifica- 
tions of the photon hopping as compared to the form as- 
sumed in equation (|5]l. A challenge here would be to pre- 
vent that significant photon population is localised in the 
fibre instead of in the cavities themselves. Nonetheless, 
even if the fibre contains appreciable excitations, interest- 
ing many-body physics is expected to happen as the fibre 
would merely play the role of an additional lattice site with 
different local terms than the other sites. 

Figure|28]shows an array of 10 micro disc cavities that 
are all coupled by one tapered optical fibre and mounted on 
an atom chip. The purpose of the atom chip is to trap atoms 
close to the cavity surfaces and, in that way, maintain a 
strong coupling regime for substantial times. 



A second class of micro-cavities that are routinely pro- 
duced in large arrays and which have very high Q-factors 
are silicon structures of either a disc 1591 or a toroidal ll57l 
shape, see figure|27] These cavities are closed silicon discs 
or toroids that typically are carried by a thin pillar in their 
centre. They are routinely produced in large arrays 11051 . 
The light is trapped in whispering gallery modes that are 
localised close to the outer surface of the structure and 
have a small mode volume II 10611 . Very efficient and highly 
tunable coupling of photons into and from those cavities is 
possible via tapered optical fibres that are placed close to 
their surface II107I . Here the evanescent field of both cav- 
ity and fibre overlap and photons can hop between both. 
The photon hopping rate can thereby be controlled very 
accurately via the separation between cavity and fibre. 

These cavities can interact with atoms in a strong cou- 
pling regime with cooperativity factors ^ ^ 50 1581 if the 



8.3. Fibre based cavities 

Strong coupling experiments with many atoms in one cav- 
ity and at the same time large cooperativity factors for each 
individual atom have been performed with fibre Fabry- 
Perot cavities i64j . These cavities are formed by two tips 
of optical fibres that are facing each other and are mounted 
on an atom chip to trap an manipulate the atoms that inter- 
act with the cavity modes. Current chips contain two such 
cavities as shown in figure |29] and cooperativity factors as 
high as ^ = 145 have been achieved l64l . 

Photons are coupled into and from the cavity directly 
via the optical fibres. Hence as for fibre coupled micro 
discs or toroids this may lead to modifications of the form 
of the photon hopping as compared to equation (|5]l. On the 
other hand, the atoms are here trapped on the atom chip 
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Figure 28 (a) Top-view optical image of a fibre taper aligned 
with and array of ten micro-disks (the remaining four micro-disks 
are out of the field of view of the image), (b) SEM image of the 
array of ten micro-disks, (c) Fibre taper transmission vs wave- 
length when the taper is aligned with the micro-disk array in (a). 
The letters match specific micro-disks in (a) with the correspond- 
ing resonances in (c). Reprinted with permission from 1591 . 



a be 




Figure 29 a, Layout of an atom chip that contains two fibre Fab- 
ryPerot cavities, b, Close-up view of the two fibre FabryPerot 
(FFP) optical cavities that are mounted on the chip. Cavity modes 
are drawn to scale in red. The BEC is produced in a magnetic trap 
and positioned in the FFPl mode ('B'). c, Geometry of the FFPl 
cavity, d, Overlay of three CCD time-of-flight (TOF) absorption 
images, showing the anisotropic expansion of a BEC having in- 
teracted for 50 ms with the cavity field. The optical fibres are 
outlined for clarity. Reprinted with permission from 1641 . 



and can thus be held in a location of strong interaction with 
the cavity modes for a long time. 



8.4. On-chip Fabry-Perot cavities 

Micro-fabricated arrays of Fabry-Perot type cavities with 
large cooperativity factors have been presented in 010911 . 
Each plano-concave resonator consists of an etched spher- 
ical mirror and a coated single-mode optical fibre, allow- 
ing open access to the intra-cavity field mode. Here, large 




Figure 30 Left: Single on-chip Fabry-Perot cavity formed by an 
etched micro-mirror and a tip of an optical fibre. Right: Two di- 
mensional array of cavity micro-mirrors. Reprinted with permis- 
sion from II 101 . 



regular arrays of spherical mirrors can be created with con- 
trollable sizes and spacings. These cavities which have a 
very small mode volume are illustrated in figure [30l The 
extremely small mode volume leads to a large single-photon 
Rabi frequency, and light is conveniently coupled into and 
out of the cavities via the fibres themselves. These cavi- 
ties to date achieve a finesse F ~ 5000, corresponding to 
Purcell factors (or cooperativities) ^ ^ 40. 

As with the other cavities that involve optical fibres 
modifications of the photon hopping as compared to equa- 
tion (|5]l are expected here as well. Current strong coupling 
experiments have also been only performed with falling 
atoms 1631 . 



8.5. Superconducting stripline resonators 

A promising experimental platform that operates in the 
microwave regime are cavity QED system formed by a 
Cooper-pair box coupled to a superconducting transmis- 
sion line resonator. A Cooper-pair box is formed by two 
superconducting islands separated by a Josephson junc- 
tion. Although it is a macroscopic device, it behaves in 
many ways like an atom as the relevant structure can be de- 
scribed by a two level quantum system. A superconducting 
transmission line cavity, on the other hand, is a quasi-one- 
dimensional coplanar waveguide resonator, formed on a 
chip with a resonance frequency in the microwave range. 
The strong-coupling regime in such a system has been achieved 
ll62l and recently two Cooper-pair boxes have been strongly 
coupled to the cavity mode II111II112I . Very large cooper- 
ativity factors have been achieved in this set-up 162111 121 
and the nonlinearity of the Jaynes-Cummings model has 
been measured 11 1311 . Although coupling between differ- 
ent cavities has not been realised, up to ten Cooper-pair 
boxes interacting in the strong coupling regime with the 
resonant cavity mode has been predicted to be achievable 
using the architecture of Ref . Ill 121 . 
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Figure 31 a, A superconducting niobium coplanar waveguide 
resonator of 24 mm length. It is coupled by a capacitor at each 
end of the resonator (see b) to an input and output feed line, b, 
The capacitive coupling to the input and output lines, c, False 
colour electron micrograph of a Cooper pair box (blue) fabri- 
cated onto the silicon substrate (green) into the gap between the 
centre conductor (top) and the ground plane (bottom) of a res- 
onator (beige) using electron beam lithography and double angle 
evaporation of aluminium. Reprinted with permission from 1621 . 



9. Applications 

Structures of coupled arrays of micro-cavities have sev- 
eral interesting potential applications. Particularly useful 
in this respect is the individual addressability of single 
cavities with optical lasers. The feasibility of single site 
addressing is based on the size of the cavities which can 
easily be larger than an optical wavelength. Furthermore 
the inter-cavity coupling could be realised by optical fibres 
which would move the cavities even further apart. 

Among the potential applications are the implementa- 
tion of tools for quantum information processing. These 
include the generation of entangled states such as cluster 
states |[27i .52,54,55il (see section lT4b . the controlled prop- 
agation of polaritonic excitations in cavity arrays II53II35I 
[36ll37ll50l . quantum gates ISTI and the generation of single 
photons in parallel 1561 . 

Perhaps most importantly, the local controllability of- 
fered by coupled cavity arrays can also give rise to novel 
physical effects and we will mention here one example of 
these. Local manipulations of an effective many body sys- 
tem could be applied to engineer a boundary between a 
Mott insulator and a superfluid region of a Bose-Hubbard 
model at unit filling. Initially both regions are decoupled 
and cooled to their respective ground states. This situation 
can for example be prepared in a cavity array that is oper- 
ated in the way as discussed in section |4] and thus hosts a 
polaritonic Bose-Hubbard Hamiltonian dTji. After switch- 
ing on a small tunnelling rate between both regions, all 
particles of the Mott region migrate to the superfluid area. 
This behaviour is shown for a one dimensional chain in 
figure |32] where the particle densities for all lattice sites 




Figure 32 The particle density across the chain as function of 
time. Sites 1 to 20 start in a Mott insulator regime and the rest in 
a superfluid regime. U = 1.0, J = 0.1, fi = 0,U = 0.2, fl = 0, 
J = 1.0 and Jj — 0.1. Reprinted with permission from 1331 . 



are shown as a function of time. The migration takes place 
whenever the difference between the chemical potentials 
of both regions is less than the maximal energy of any 
eigenmode of the superfluid Ii33il . 



10. Summary 

In this review, we have tried to give an overview of the cur- 
rent status of approaches to engineer effective many body 
dynamics in arrays of cavities that are coupled via photon 
tunnelling. The schemes that have been developed so far 
can be divided into two main categories, effective Bose- 
Hubbard type models and effective spin models. 

In Bose-Hubbard type models, quantum phase transi- 
tions between superfluid phases of polaritons or photons 
(coherent light) and Mott insulators can be studied. Here, 
the experimentally challenging and hence physically new 
and interesting part is to create a Mott insulator in which 
polaritons or photons are by their effective mutual repul- 
sion forced to stay in one lattice site (cavity) and not able 
to move. This situation would in correspond to a crystal 
formed by light or "frozen light". 

A key issue on the experimental side is to generate 
a Kerr nonlinearity or effective interaction between po- 
laritons/photons that is significantly larger than the decay 
mechanisms present in the cavity array. 

Existing approaches for effective spin models mainly 
employ two metastable levels of an atom trapped in each 
cavity to represent the two spin polarisations f and j. Ra- 
man transitions mediated by external driving lasers that 
give rise to interactions between neighbouring cavities via 
the exchange of virtual photons are then used to generate 
effective spin-spin interactions. 

As for Bose Hubbard type models a main experimen- 
tal challenge is to make atom-photon interactions signif- 
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icantly stronger than the decay mechanisms in the cavity 
array. 

Effective many body system in a coupled cavity arrays 
are of great interest from various points of view. First of all 
they open up a way to study condensed matter physics with 
light and hence to engineer new quantum states of light. 
Furthermore they offer a great versatility for applications 
as quantum simulators due to the possibilities to manip- 
ulate them via external driving lasers. Most importantly 
they would, by construction, allow for control, manipu- 
lations and measurements of individual lattice sites. The 
situation is somehow complementary to cold atoms in op- 
tical lattices, where very large lattices can be created and 
experiments provide good measurement access to global 
quantities such as the phase coherence of the superfluid. 

A polaritonic many body system in a cavity array on 
the other hand will for the foreseeable future only com- 
prise much smaller lattices, but provide good access for 
the measurement of local quantities. In contrast to optical 
lattices it would thus allow to directly verify the integer oc- 
cupation number of each lattice site in the Mott insulator 
regime. 

The local addressability of coupled cavity arrays gives 
rise to a number of novel application. These contain quan- 
tum information processing application but also novel phys- 
ical situations that give rise to features such as particle mi- 
gration. 
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